Life of the nodal quasiparticles in Bi-2212 as seen by ARPES 
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While the pronounced doping dependence of the quasiparticle spectral weight in the antinodal 
region of the superconducting cuprates, as seen by ARPES, unambiguously points to the magnetic 
origin of the strong electron-boson coupling there, the nature of the electron scattering in the nodal 
direction remained unclear. Here we present a short review of our recent detailed investigations of 
the nodal direction of Bi-2212. Our findings prove the existence of well defined quasiparticles even 
in the pseudogap state and show that the essential part of the quasiparticle scattering rate, which 
appears on top of Auger-like electron-electron interaction, also implies a magnetic origin. 

PACS numbers: 74.25.Jb, 74.72. Hs, 79.60.-i, 71. 15. Mb 



I. ARPES VIEW 

Angle-resolved photoemission spectroscopy (ARPES) 
provides a direct view on the density of low energy 
electronic excited states in solids — the 2D detector of 
the electron analysers used in modern ARPES is just a 
window into momentum-energy space of 2D compounds. 
A snapshot through this window stores the quasiparti- 
cle spectral weight in the momentum-energy co-ordinates 
H, H, ^ 1^. Being essentially two-dimensional, the su- 
perconducting cuprates are a perfect example of the 
"arpesable" compounds 0|. All the interactions of the 
electrons which are responsible for their unusual normal 
and superconducting properties are encapsulated in such 
snapshots, and success in understanding of the nature 
of electronic interactions in the cuprates depends, in the 
first place, on how clear the ARPES window is. Then, 
the experimental experience (namely, how many differ- 
ent snapshots have been taken and made out) comes into 
play. But taking into account a number of parameters 
(e.g. temperature and doping) which cause a redistri- 
bution of the quasiparticle spectral weight, the detailed 
exploration of the momentum-energy space, even for one 
compound, will take ages of experimental work. 

Leaving such a global task for the nearest future, one 
can focus on the two cuts in the Brillouin zone (BZ): 
nodal and antinodal directions (see Fig.^a). These re- 
gions represent an inherent anisotropy of the electronic 
interactions in the cuprates which appear in anisotropy 
of the superconducting gap 0, pseudo-gap and cou- 
pling strength ^ (or scattering in general). While the 
pronounced doping dependence of the quasiparticle spec- 
tral weight in the antinodal region of the BZ unambigu- 
ously points out to the magnetic origin of the strong 
electron-boson coupling seen by ARPES H ES 113 (see 
Fig-fflc-f), the nature of the electron scattering in the 
nodal direction remaines unclear. Here we report the 
results of a detailed investigation of the nodal direction 
of Bi-2212 in a wide range of doping, temperature and 
excitation energy. We have found that although the elec- 
tronic band structure along the nodal direction remains 
complex due to non- vanishing bilayer splitting, the quasi- 



particle spectral weight distribution from each split band 
can be self-consistently described within the quasiparticle 
self-energy approach. The scattering rate, on the other 
hand, can be considered as a sum of two main channels: 
the doping independent channel can be well understood 
in terms of the conventional Fermi liquid model, while the 
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FIG. 1: Electronic band structure of an overdoped Bi-2212 
ii]: Fermi surfaces (a) and the "XMY" cut (b). ARPES 
snapshots taken along the XMY direction of the BZ fiot for 
the overdoped (c, d) and underdoped (e, f) samples above (c, 
e) and below (d, f) T^. 
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additional doping dependent channel implies a magnetic 
origin. 



II. EXPERIMENTAL CORNERSTONES 

The experimental details can be found elsewhere 0, 
but here we highlight the cornerstones which are 
peculiar to our experiments. They are: the precise cryo 
manipulator, the wide excitation energy range, and the 
superstructure-free samples. 

The precise cryo manipulator operates in the controlled 
temperature range from 20 K to 400 K and allows us 
to translate the sample in three perpendicular directions 
and rotate it around three perpendicular axes in steps of 
0.1°, that secures the precise positioning and easy motion 
of the ARPES window in the momentum space. 

The photons of different energy and polarisation have 
appeared to be an extremely useful tool to selectively 
excite the electrons from different bands 0,111. As a 
light source we use a He discharge lamp, linearly polar- 
ized light from a high resolution beamline (U125/1-PGM 
at BESSY) with a wide excitation energy range (17-600 
eV), or circularly polarized light (4.2R beamline "Circu- 
lar Polarization" at ELETTRA). 

The well known problem for ARPES on Bi-2212 is 
a "5x1" superstructure which produces a number of 
diffraction replicas along the BZ diagonal (FY direction) 
[T^ITsl l. This highly complicates the analysis of the spec- 
tra taken from certain areas of the BZ (e.g., the antinodal 
region 12, 16J) or, if the main features and replicas are 
spatially separated (e.g. along FY direction), reduces the 
photocurrent intensity. Another important point about 
Bi-2212 samples, which becomes crucial now is 
that one needs to know precisely their doping level. We 
studied both the led-doped superstructure-free Bi(Pb)- 
2212 and reach of replicas pure Bi-2212. However, as the 
main line, we use the well characterized superstructure- 
free samples of a wide doping range (x = 0.11 - 0.22), 
for which the charge carrier densities have been derived 
from the measured Fermi-surface area (and appeared to 
be consistent with their Tc) and the tight-binding 
parameters have been determined 6] . The parameters of 
the band structure of the pure Bi-2212 samples which we 
have measured are also in agreement with their doping 
level estimated from the Tc measurements. 



III. COMPLEX STRUCTURE 

A distinguishing feature of modern ARPES is the abil- 
ity to resolve the bilayer splitting (BS) of the CuO con- 
duction band in the bilayer cuprates. For the first time 
such a splitting has been observed for overdoped Bi-2212 
[20I I21I and then also for optimally doped and under- 
doped sam ples [l9l l2^ (clearly resolved below 0, l2^ 
and above llfl the superconducting transition). It has 
been found ^20l I2H that the observed splitting can be 
approximated by a momentum dependence: t±{coskx — 




FIG. 2: Zoomed ARPES snapshots taken along the nodal 
('TY") direction of the BZ for Bi-2212 UD80 (a,b) at 27 eV 
and 17.5 eV excitation energy respectively |25ll . The experi- 
mental data for the same sample presented in form of energy 
distribution curves (EDCs) (c) and momentum distribution 
curves (MDCs) (d). EDCs are taken in the momentum range 
from kp- 0.025 A"^ (top) to kp + 0.015 A"^ (bottom), where 
kp is an average between antibonding, kp, and bonding, kp, 
Fermi level crossings; red EDCs roughly correspond to kp and 
kp. MDCs are taken in the energy range from 3 meV (top) 
to -27 meV (bottom); Ep-MDC shown in red. (e) MDCs in- 
tegrated in energy (along the experimental dispersion) about 
10 meV from Ep: (1) - Bi-2212 UD80, 1st BZ, 17.5 eV; (2) - 
Bi-2212 OP89, 2nd BZ, 20 eV for the bold curve and 18 eV for 
the dotted curve; (3) - Bi(Pb)-2212 OD73, 1st BZ, 17.5 eV, 
dashed curve results from the same MDCs but normalized to 
highest binding energy; (4) - Bi-2201, 1st BZ, 17.5 eV; the 
red curves, when shown, represent fitting results. 



cosfcy)^/2, which is expected for an inter-plane hopping 
between two CuO layers (where t± describes the inter- 
layer hopping mainly mediated via Cu4s orbitals). The 
splitting along the node is expected to be not zero but 
much less than the maximum splitting at the saddle- 
point jl^. In order to answer the question whether the 
splitting really vanishes in nodal direction, we performed 
precise measurements in the low excitation energy range 
(17-22 eV) [23. The total energy resolution was set to 
10 meV, the angular resolution of the analyser was 0.15°. 
To conclude on the existence of the nodal BS one should 
ensure taking the spectra from exactly the nodal direc- 
tion. We determined the nodal direction measuring the 
Fermi surface (FS) maps with 0.5° step in azimuth an- 
gle (for the details about the experimental setup see 5]). 



3 



The spectra which we quahfy as nodal and discuss below 
are taken from the FS cuts with the smallest kp which 
also turned out to have the steepest dispersion and the 
smallest leading edge gap 0. 

Fig. 12 represents the experimental evidence for the 
nodal splitting. Panels a and b show ARPES snapshots 
taken along the nodal (FY) direction of the BZ for Bi- 
2212 UD80 sample at 27 eV and 17.5 eV excitation energy 
respectively (note that the momentum-energy window is 
much smaller here than on standard snapshots like in 
Fig.^l. While at 27 eV one can see only one band cross- 
ing the Fermi-level, two bands are clearly visible at 17.5 
eV. One can also notice the presence of two bands in the 
energy distribution curves (EDCs), see panel c, extracted 
from the same dataset and, more explicitly, on the mo- 
mentum distribution curves (MDCs) presented in panel 
d. In panel e (curve 1), in order to improve statistics, we 
integrate the MDCs along the experimental (renormal- 
ized) dispersion in the energy range 10-20 meV around 
Ep, where the MDC width does not vary dramatically. 
Panels a-d show one example but we observe the same 
effect on a number of samples of different doping level, 
with and without 5x1 superstructure . 

The dependence of matrix elements on excitation en- 
ergy for the nodal point in the 1st BZ exhibits a local 
maximum at about 17.5 eV for both the total intensity 
from bilayer split band and the intensity from the bond- 
ing band compared to its antibonding counterpart. For 
the nodal point in the 2nd BZ the dependence on matrix 
elements is different and the bonding band is the most 
pronounced for hiy = 20-21 eV excitation energy: the 
MDCs 2 in Fig. show how the bonding band peak 
appears when going from 18 eV (dotted curve) to 20 
eV (bold curve). The observed excitation energy depen- 
dence of the effect is in accord with recent calculations of 
ARPES matrix elements |2^ which show that at low en- 
ergy range the emissions are dominated (peaked at about 
18 eV) by excitation from just the O sites. 

To extract precise BS values we fit the integrated 
MDCs to a superposition of two independent Lorentzians 
(an example of a fitting curve is shown in Fig. [2|3) . For 
the presented dataset the splitting in momentum Afc — 
0.012(1) which corresponds to 48(4) meV bare band 
splitting (for bare Fermi velocity vp — 4.0 eVA 0,|23|) or 
23 meV splitting of the renormalizcd band (renormalized 
Fermi velocity Vp = 2.0 eVA). For other bilayer samples 
the values are similar and in a good agreement with the 
LDA band structure calculations [l^ . This enables us to 
assign the splitting predominantly to vertical inter-plane 
hopping between 02p„ orbitals and to conclude on the 
lack of any electronic confinement to single planes within 
a bilayer in Bi-2212 due to strong correlations. We note 
that such a careful comparison between theoretical and 
experimental values cannot be done anywhere in the Bril- 
louin zone except the node due to gap opening (along the 
nodal direction, analysing MDCs, we determine the dif- 
ference in kp for these two bands while in other fc- regions 
one can mainly rely on EDC analysis in which the rela- 
tion of EDC peak position with the position of the band 



is highly model dependent, e.g. see 0|). 



IV. SIMPLE PHYSICS 

Apparently, the nodal splitting found in Bi-2212, if not 
taken into account, can intricately complicate the nodal 
spectra. Even if unresolved, it should influence the quan- 
tities derived from these spectra, e.g. the renormalized 
dispersion or scattering rate. It is observed as a depen- 
dence of these quantities on excitation energy. In the 
following we try to eliminate the influence of the split- 
ting choosing an appropriate excitation energy. We focus 
on the experimental dataset taken at 27 eV in order to 
find out whether the quasiparticle spectral weight from 
only one band can be described by "simple physics", i.e. 
in terms of quasiparticle self-energy . 

Fig.Oillustrates the basics of the nodal spectra analy- 
sis within the self-energy approach. The black solid line 
represents a non-interacting case when the spectral func- 
tion is a delta function with the pole u) — ei}^^) — 0. 
£(k) is called "bare dispersion". A simple electronic in- 
teraction leads to shifting and broadening of the non- 
interacting spectral function and the resulting picture 
is essentially that which is measured by ARPES: the 
blurred region in Fig. illustrates the distribution of the 
quasiparticle spectral weight. Such an interaction can 
be described by introducing a quasiparticle self-energy 
Y, — Yi' + iT," , an analytical function the real and imag- 
inary parts of which are related by the Kramers-Kronig 
(KK) transformation . Neglecting the momentum de- 
pendence of the self-energy, the MDC at certain ui ex- 
hibits a Lorentzian lineshape 4] with the maximum at 
kmii-o) determined by S'(a;) — lo — e{km), which is illus- 
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FIG. 3: Bare band dispersion (solid line) and renormalized 
dispersion (red points) on top of the spectral weight of inter- 
acting electrons. Though intended to be general, this sketch 
represents the nodal direction of an underdoped Bi-2212 |2'^. 



4 



trated in Fig. O by the double headed arrow (red squared 
symbols show the renormalized dispersion). In the re- 
gion where the bare dispersion can be considered as lin- 
ear {e = VFk), the MDC width W (the half width at half 
maximum) is proportional to S": S"(w) = —vfW{uj). 
Thus, the determination of both the real and imaginary 
parts of the self-energy requires the knowledge of the bare 
dispersion. The KK transformation, giving an additional 
equation to relate these functions, opens the way to ex- 
tract all desired quantities from the experiment. 

For example, one can express the coupling strength 
A = —{dYi'{uj)/(hj)^=Q = vf/vr — 1, where vn is the 
renormalized Fermi velocity, as |27l | 



A 



-DE". (1) 



Using the above definition of the D operator, Vp^ — 



DW, or 1 + A = where 
Z =1- vrBW 



(2) 



is the coherence factor (0 < Z < 1). 

In case W{uj) decays to zero or saturates on the scale 
covered by experiment, as it is expected for the scattering 
by phonons [30j, the parameters vp, X oi Z can be eas- 
ily determined from the experimental values of vr and 
DW. In cuprates, however, the MDC width W along 
the nodal direction does not decrease or even saturate in 
the whole experimentally accessible energy region (up to 
= 0.5 eV) and one can make only a rough estimation 
expanding T>W = T>'^'"Wexp + T^Z',„Wmod, where Wexp is 
the experimentally determined function of w, and Wmod 
is a model function which depends on both the high en- 
ergy cut-off, \ujc\ > \(^m\-, above which S"(ti;) starts to 
decrease or saturate and a model for these high energy 
tails. For a simple estimation one can take Wexp — ol'^'^ 
and Wrr,.nd = ckwi which gives Tif^W^xp = D,T W, 



mod: 



demonstrating that the contribution of the tails can be 
essential. 

Fortunately, we have found that with some assump- 
tions and considering a sufRciently wide energy range 
(up to 300 meV) one can precisely determine the bare 
dispersion while the cut-off energy for E"(w) and its tails 
remain undefined. In this procedure we use the parabolic 
bare dispersion which is a good approximation of the dis- 
persion derived from the tight-binding fit of the Fermi 
surface 0. This assumption brings some corrections to 
the self-energy parts derived from the experimental data: 



I]'(c.) 



VF 



4) 



T."(uj) = -^VFW(uj)^kl^{uj) - W^{uj). 
The KK transform completes the system: 
I]'(w) = KKS"(cj). 



(3) 
(4) 

(5) 



Eqs. (0) and ||SJ| give two independent ways to calcu- 
late E'(w), and this is a core of the procedure which is 
described in details in 1271. 
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FIG. 4: Real and imaginary parts of the self-energy extracted 
from the experiment with the described procedure. A com- 
plete coincidence between the corresponding parts of the self- 
energy calculated from the two different experimental func- 
tions, the MDC dispersion and MDC width, demonstrates 
the full self-consistency of the ARPES data treated within 
the self-energy approach |2^. 



Fig.^shows the results for an underdoped Bi(Pb)-2212 
(Tc = 77 K) at 130 K, i.e. in the pseudo-gap state. S'(w) 
calculated in two ways have appeared to be identical for 
the certain vf = 3.82±0.17 eVA. This gives the following 
interaction parameters: A = 0.87 ±0.12, Z = 0.54 ±0.03. 
The extracted bare band dispersion is in good agreement 
with the band structure calculations j2|| and allows one 
to quantify the self-energy of the electronic excitations 
on the real energy scale. 

The demonstrated self-consistency can be considered 
as a validity criterion for photoemission spectra, since it 
weeds out not only a "complex structure" like splitting 
or admixture of other bands but also artificial effects like 
inhomogenity of the detector, etc. Considering the in- 
timate relation between S' and S" one can distinguish 
two energy scales in Fig.^ One, at about 200 meV, cor- 
responds to the maximum of S'(w) and is related to the 
cut-off energy for ^"(0;). Another, at 70 meV, is a famous 
" kink" , it develops as a sharp bend between two linear 
segments of S'(a;) (a peak in d^Y^' (uj)/ dixj'^Yjs related to 
the " drop" in the scattering rate [S^, 113, 13 , and com- 
monly explained as an interaction with a bosonic mode 

13 Q . It is important to distinguish these two scales 
in order to find out the nature of the coupling boson. 

As far as the kink on the dispersion appears as just a 
sharpening of a bend of the same sign in the experimen- 
tal dispersion which is present at every temperature and 
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FIG. 5: (a) Temperature dependence of the scattering rate for the nodal quasiparticles in optimally doped Bi(Pb)-2212. The 
gray solid line represents a contribution from the usual Auger decay (Fermi liquid parabola) 28] obtained by fitting the data 
for highly overdoped sample (OD69) at 130 K (see inset), (b) Result of a subtraction of the Fermi liquid parabola for each 
temperature in terms of the imaginary part of the self-energy (the FWHM/2 is multiplied to the bare Fermi velocity vf ~ 
4 eVA). (c) Strengthening of the scattering mode with underdoping. Comparison of the imaginary part of the self-energy of 
nodal quasiparticles in Bi(Pb)-2212 underdoped (Tc — 76 K) and overdoped (Tc = 73 K) samples at 25 K. The shaded areas 
represent the contributions from the magnetic scattering obtained by subtraction of the FL parabola |34|]. 



doping 1^ , we focus on the " scattering rate kink" 
|34| which is much more convenient in this sense because 
it develops on top of the strong normal state scattering 
of the opposite curvature. 

Studying a number of samples of different doping level 
at different temperature we have found that the scatter- 
ing rate kink makes it possible to distinguish between 
the different scattering channels ■ We argue that the 
main contribution to the scattering can be well under- 
stood in terms of the conventional Fermi liquid model 
(FL) |28| while the additional doping dependent contri- 
bution apparently has a magnetic origin. 

Fig.|3^ shows the scattering rate (in momentum units) 
as a function of frequency for optimally doped Bi(Pb)- 
2212 OP89 for different temperatures. A sharp kink seen 
in T,"{uj) at 0.1 eV (indicated by the arrow) at 40 K 
(below Tc = 88 K) gradually vanishes with increasing 
temperature. Another important finding is that the high 
binding energy tail of S"(w) shifts upwards with tem- 
perature similar to the £"(0) value. This sh ift, being in 
agreement with optical conductivity resul ts l35l| . contra- 
dicts, in fact, the marginal FL scenario p3a |. according 
to which E"(cj,T) oc max(|a;|,r). Such a shift of the 
whole curve is expected within the FL model when the 
scattering rate is determined by an Auger-like decay (the 
process where the hole decays into two holes and one elec- 
tron H) that gives E" oc uj'^ + {ttT)'^ . The FL behaviour 
is generally expected for overdoped samples and in 
Fig. [SJi we add the FL parabola (solid line) which per- 
fectly fits the scattering rate for an OD69 sample above 



Tc in the whole binding energy range. This parabola 
evidently describes the main contribution to E" at any 
temperature. The additional contribution, which is seen 
as a hump on top of the FL parabola, must originate 
from an additional interaction which can be responsible 
for the unusual properties of the cuprates. In Fig. I^b 
we evaluate this interaction subtracting the FL parabola 
for each temperature and setting the resulting offsets to 
zero. 

Therefore, this additional contribution decreases with 
increasing temperature, and we have found that it van- 
ishes above Tc for the overdoped samples, but persists 
at higher temperatures, presumably up to T* for opti- 
mally doped and underdoped samples [s^l- In Fig. ]^ 
we compare the absolute values of E"(w) for underdoped 
(UD76) and overdoped (OD73) Bi(Pb)-2212 at T = 25 
K. The room temperature scattering rates for these two 
samples coincide within the experimental error bars. It 
is seen that at low temperature the underdoped sam- 
ple exhibits a much higher scattering rate with a more 
pronounced kink that has a tendency to disappear com- 
pletely at higher doping levels [s^l- The differences be- 
tween these data and the FL parabola (solid line, the 
same as in Fig. la) demonstrate that the additional scat- 
tering channel of the nodal quasiparticles is highly doping 
dependent which is difficult to reconcile with the phonon 
scenario Ull , leaving space for magnetic excitations 
as the only bosons responsible for this additional channel 

[iilii- 

Although our findings support the magnetic nature of 
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the doping dependent channel in the scattering rate, its 
origin is still to be understood. There are two main sus- 
pects for the scattering kink and dispersion kink (in terms 
of a sharp feature that was discussed above): the mag- 
netic resonance and gapped spin-fluctuation continuum 
[4l| . As far as the magnetic resonance at (tt, tt) is be- 
lieved to be sharp in energy and momentum it seems 
unlikely to connect the nodal region to another part of 
the FS by such a scattering. This increases the possibil- 
ity to describe the kinks by the gapped spin-fluctuation 
continuum alone and it has been recently shown that 
the gapped magnetic spectrum can well describe the dis- 
persion kink feature. Nevertheless, we believe that the 
magnetic resonance cannot be ruled out because of the 
presence of the so called "shadow band" ^3 which, as 
we have recently shown is not a diffraction replica, 
but a real band of the CuO layers, and consequently pro- 
vides an easy way for the quasiparticle from the main 
band to scatter. In order to find out which mechanism is 
dominant, more accurate and systematic investigations 
of both kink features are needed. 



V. CONCLUSIONS 

We performed a detailed investigation of the nodal di- 
rection of Bi-2212 in wide ranges of doping, tempera- 



ture and excitation energy. Even along this direction, 
the electronic band has appeared to be split, but choos- 
ing an appropriate excitation energy we have succeeded 
to single out the quasiparticle spectral weight distribu- 
tion from only one bilayer split band and show that it 
can be self-consistently described within the quasiparti- 
cle self-energy approach. Focusing on the scattering rate 
as a function of binding energy, temperature and doping 
we have distinguished two main channels in the electron 
scattering. While the main doping independent channel 
can be well understood in terms of the conventional Fermi 
liquid model, the additional doping dependent channel 
implies a magnetic origin. Our findings prove the exis- 
tence of well defined quasiparticles even for the under- 
doped Bi-2212 (Tc ~ 77 K) in the pseudogap state. 
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